Domain instability during precessional magnetization reversal 
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Spin wave equations in the non-equilibrium precessing state of a ferromagnetic system are found. 
They show a spin-wave instability towards growing domains of stable magnetization. Precession of 
the uniform magnetization mode is described by the Landau Lifshitz equation with the exponentially 
growing in time effective Gilbert dissipation constant that could have both signs. On the developed 
stages of the domain instability a non-stationary picture of domain chaos is observed. 
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Recently fast pulses of strong magnetic field have 
brought a possibility to study experimentally a preces- 
sional magnetization reversal 0, 12, @] ■ The magnetic 
field of such a pulse is spatially uniform on the lengths 
that are typical for a ferromagnetic system: the domain 
wall thickness, the grain size or the width of thin ferro- 
magnetic films etc. The amplitude of the magnetic field 
in the pulse is such that the Zeeman energy exceeds the 
anisotropy energy that is present in the ferromagnetic 
system. The integral of the magnetic field over the du- 
ration of the pulse could be varied to result in a uniform 
rotation of spins by angles ranging from few degrees to up 
to 2ir pj. There are different possibilities how to apply 
the magnetic field vector with respect to the anisotropy 
axis of a ferromagnetic system Q ■ Also little change oc- 
curs to the ferromagnetic order of rotating spins during 
the pulse. All these special properties of the fast pulses 
allows one to prepare in effect a non-equilibrium uniform 
ferromagnetic state with the magnetization vector point- 
ing in an arbitrary direction. An alternative method to 
reverse magnetization is an application of longer pulses 
of weaker magnetic field where magnetization precession 
results from the combined action of the magnetic field 
and the anisotropy during the pulse El El El ■ 

In general after the pulse ceases a precession of uniform 
magnetization will start due to the anisotropy with the 
magnetization unit vector sweeping large angles. There is 
an emerging experimental evidence that such a precession 
follows the Landau Lifshitz equation with the Gilbert dis- 
sipation constant exceeding by tens times the Gilbert dis- 
sipation constant in the same ferromagnetic system un- 
der the conditions of the Ferro-Magnetic Resonance 
where small angles of precession are excited. There are 
two special non-equilibrium ferromagnetic states. First, 
in some systems the dissipation is especially weak and 
the large angle precession is periodic in time. Second, 
the magnetization unit vector after the pulse is rotated 
close to the maximum or saddle point of the anisotropy 
energy. Obviously, in the last case a fundamental insta- 
bility towards growing domains of stable magnetization 
arises with the domain sizes being of the order of the 
domain wall thickness. The initial stages of this domain 
instability can be described in terms of the growing spin 



waves. In the first case in general a spin wave instability 
also develops like the two-magnon Suhl instability that 
has been extensively studied in connection with the para- 
metric excitation of spin waves in a ferromagnet in equi- 
librium 01 • Recently such instability has been shown to 
be responsible for a large effective damping during peri- 
odic precession in the area of angles around the stability 
point [13. 

In this letter I find a system of spin wave equations in 
the rotating frame that follows the uniform mode preces- 
sion of the non-equilibrium state of ferromagnetic sys- 
tem. Solution of these equations under certain condi- 
tions yields a spin- wave instability. The equation for uni- 
form mode precession coupled to the dynamics of grow- 
ing instable spin waves is the Landau-Lifshitz equations 
with effective Gilbert constant that exponentially grows 
in time too. On the developed stages of instability ab- 
initio numerical solution of the Larmor precession in a 
cluster of exchange coupled classical spins shows a non- 
stationary picture of domain chaos. 

Consider a ferromagnetic particle in a precessional 
state that is described by the uniform time-dependent 
ferromagnetic order parameter unit vector n(t). Spin 
waves are also present in this particle either due to the 
thermal excitation at temperature T before arrival of the 
pulse or due to the spin wave instabilities developing after 
the pulse. The local ferromagnetic order parameter unit 
vector M(t,r) = n(t)y/l -4> a 4> a + (t> a e a {t) (M 2 = 1) is 
parameterized by spin wave two-component field 4> a (t, r) 
with a — 1,2 and by a rotating frame of three mutu- 
ally orthogonal vectors n, e a at any moment of time. 
A specific choice of two vectors e a (t) is dictated by the 
vector n(t) apart from the arbitrary time-dependent ro- 
tation about the axis n. This rotation represents U(l) 
gauge invariance in the description of spin waves in a fer- 
romagnet. It is convenient to relate the rotating frame 
to the curve w(t) in the spin space. The vector n(t) 
is the unit tangent vector to this curve: dw/dt = n(t), 
the vector e%(t) being the unit principal normal vector 
to this curve and ea(t) being the unit principal binor- 
mal vector to this curve. This choice imposes a gauge 
fixing condition known as the Frenet-Serret equations: 
dn/dt = (xn, de a /dt = £ x e a , where the rotating frame 
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is precessing as a solid structure about the Darboux vec- 
tor = —nri+ ke2, where k(t) is the local curvature 
of the curve and n(t) is the local twist of the curve asso- 
ciated with the connection of the U(l) gauge. 

We consider low temperatures with only long wave- 
length spin waves being excited. Because the pulse rep- 
resents a fast coherent rotation new spin waves with short 
wavelength are not exited during the pulse. Therefore the 
induced spin precession afterwards is classical long wave- 
length dynamics. The exchange part of the Hamiltonian 
of the ferromagnetic particle reads 
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r, 
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where J is the exchange constant and d is the dimen- 
sion of a ferromagnetic system. The Goldstone theorem 
forbids contribution to the exchange Hamiltonian from 
the uniform mode n. Relativistic spin-orbit and dipole- 
dipole interactions induce the anisotropy energy that for 
simplicity we write in the local approximation: 



H an [M] 



E ( M(t, r) ) d d f, 



(2) 



where the anisotropy density E(m) is an arbitrary func- 
tion of m (here and below m is an arbitrary unit vec- 
tor distinct from the special time dependent n(t)) and 
does not depend on the coordinate r. One-ion magnetic 
anisotropy is an even function E(m) = E(—m) and is 
symmetric under the crystallographic group of transfor- 
mations Tf. E(Tifn) = E(fh). But a uniform magnetic 
field could be included also. The choice of anisotropy as 
a uniform function of local ferromagnetic order Eq. Q is 
very general and applies to many experiments. It holds 
for particles of special shapes like the ellipsoid where the 
dipole-dipole interaction is described by effective uniform 
demagnetizing field. Also on the interfaces Neel predicted 
a quite general phenomenon of the interface anisotropy 
energy due to the break down of the translational sym- 
metry in the direction perpendicular to the interface. In 
this case Eq. @ holds provided the spin wave modes with 
wave vectors perpendicular to the interface are not ex- 
cited as for ferromagnetic thin films. 

The dynamics of a ferromagnetic particle is described 
by the classical Larmor precession equation for local fer- 
romagnetic order M(t, r): 



d t M = 7 M x 
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where 7 = e/mc is the gyromagnetic ratio and H = 
H e x[4>] + H an [M] is the Hamiltonian. The Larmor equa- 
tion conserves the total energy dH/dt — 0. We substitute 
into Eq.0 5H ex /SM = e a 5H ex /5(j) a and 8H an /8M = 
—F(M), where the force F(fh) = —dE/drh is a vector 
field on a sphere fa 2 = 1. In the Frenet-Serret gauge for 



the rotating frame: e*2 • F = k(t) and e\ -F = 0. The total 
energy is the sum of the uniform mode energy E(n) and 
the spin wave energy £ sw (n, [(f)]): H = E(n)+£ sw (n, [(f)]). 

We project the Larmor equation Eq.J2J onto an equa- 
tion representing the uniform mode Eq.Q and onto an 
equation representing the remaining non-uniform modes 
Eq.(|SJ|. The component of both uniform and non- uniform 
equations that is parallel to n(t) is a linear combination 
of the other equations and therefore is redundant. We 
assume that the spin wave modes are weakly excited and 
expand the Larmor equation up to the second power of 
spin wave amplitudes. We find the equation for the uni- 
form mode: 
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Here £ sw (n, [<fi]) is the spin wave energy in the harmonic 
approximation: 
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\ J (j(d^ a ) 2 + <f> a C ab cf> b -n-F^ a ^d d r, 



where the anisotropy energy mass tensor is: 
r< r->\ a d 



(5) 
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The gauge dependent component of the force n ■ F and 
C\\ are related to the precession twist n{t). In the limit 
of zero spin wave amplitudes the conservation of energy 
E(n) in the rotating Frenet-Serret frame gives: n ■ F — 
Cn (ft) — n(t). C22 (n) — C11 (ft) and C\% (ft) are both 
gauge independent. Obviously, the uniform mode Eq.Q 
has an integral: E(ft) +£ sw (n, [<f>]) = const, which is the 
total energy conservation. 

Note that Eq. Q has the form of the Landau Lifshitz 
equation with the Gilbert dissipation constant 



G(t) 
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k 2 (t) dt 



£ sw (n, [(/)}) 
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The sign of the Gilbert constant can be established in 
the vicinity of the maximum or minimum point of the 
anisotropy energy E(n) under the condition of growing 
amplitudes of spin waves. In the vicinity of the mini- 
mum point both the curvature k(t) > and the spin- wave 
mass tensor C a b(t) approach zero and are proportional to 
the angle of the deviation from the minimum point. On 
the contrary the twist n(t) approaches some finite neg- 
ative value. In Eq.JSJ we retain only the leading large 
term £ sw = -«(*) J c/) a (t,r)(j> a (t,r)d d f. For growing spin 
wave amplitudes and negative k: we find d£ sw /dt > 0. 
The Gilbert constant is positive and due to the condition 
E(ft) +£ sw = const: dE(n)/dt < 0. Therefore the vector 
ft approaches the minimum point during the precession. 
Superficially this could be seen as the dissipation pro- 
cess although the non-uniform Larmor precession in the 
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ferromagnetic system conserves energy and is fully time- 
reversible. The same analysis shows that near the max- 
imum point the twist n{i) approaches some finite posi- 
tive value and therefore the Gilbert constant is negative, 
dE(ri) / dt > and the vector ri approaches the maximum 
point during the precession that is accompanied by the 
spin-wave instability. 

We find also the equation for the non-uniform mode: 

eabdt^ = -Jdfya + C ah {n)ct> b +{k-h-F) 0°, (8) 

or in the principal normal and binormal coordinates: 

dtfa = - t /a^ 2 + (C 22 (n)-C , 11 (n))0 2 + C 21 (n)0 1 
-d t <t>2 = -Jdfa + C 12 {n)<t> 2 . (9) 

This system is linear and Hamiltonian. The coefficient 
Ci 2 (t?) can be related to the curvature k(t). We note 
that kei = [dE/dri x ri] = dri/dt and ke 2 = dE/dri — 
ft(ri • dE/dri), and therefore 

„ , ^ 1 f dri ti \ d , , . , , 

^^l**)^^^' (10) 

In the limit of infinitely small amplitudes of all spin 
waves the anisotropy energy E(ri) is approximately con- 
served and the uniform mode precesses periodically in 
time ri(t) = ri(t + T), where T is the period of preces- 
sion. Solution of the Eqs.© for a spin wave mode with 
momentum q can be formally written as a monodromy 
matrix: 

$(t + T)=M(t,T,q)$(t) (11) 

where <p(t) = (</>i(t), (f> 2 (t)). Monodromy matrix has unit 
determinant det M(t,T,fl = 1. Let a = TrM/2. If 
\a\ < 1 then the eigenvalues of M. are two complex num- 
bers Ai,a = exp(±i arccos(a)) otherwise (\a\ > 1) they 
are real: Ai,2 = exp(±arccosh(|a|)). In the later case 
spin wave mode is instable and its amplitude grows ex- 
ponentially with the increment v{q) — log | Ai | . Two 
examples have been studied numerically: a) E{m) = 
A^m 2 , + ^yiriy + A z to 2 and b) E(fri) = m\ + rriy + m\. 
In the case a) the uniform precession: ri(t), e*i(i), e 2 (i) 
could be found analytically in terms of the elliptic func- 
tions. Numerical solution in both cases a) and b) shows 
that a(A) is oscillating (but not periodic) function of 
spectral parameter A — Jq 2 . Note, that for C a b = 0, 
a(A) = cos(AT). There are infinitely many intervals of 
instability (zones |a(A)| > 1) in both cases, they ap- 
pear periodically and for A — > oo they shrink. We pro- 
pose that for any initial condition precession in arbitrary 
E(rri) there exists unstable spin-wave modes in momen- 
tum space. 

In the extremum point the uniform mode ri is constant 
and the curvature tensor C a b does not depend on time. In 



the adiabatic limit of very slow precession in the vicinity 
of extremum point the spin wave dispersion reads: 

Oj(q) = ^Jq 2 (Jq 2 + C 22 - C n ) - C 2 2 . (12) 

The spin wave instability develops if either a) Ci 2 ^ or 
b) C 22 — C\\ < 0. According to the condition b) systems 
with an easy-axis anisotropy always instable whereas 
easy-plane systems could be stable provided Ci 2 = that 
according to the Eq. (|10f> requires a precession along a cir- 
cle with the constant curvature. Thus a paradoxical situ- 
ation arises. Let excite a however small angle precession 
of easy-axis ferromagnet around the equilibrium point at 
such low temperature that amplitudes of thermal spin 
waves are much less than the angle of precession. Ac- 
cording to Eq. I|12l) and according to our numerical simu- 
lation of very long (2048 sites) one dimensional ferromag- 
net chain a spin-wave instability will develop. But spin 
wave dispersion in the laboratory (non-rotating) frame is 
ui(q) = Jq 2 + A, where positive A is the strength of easy- 
axis anisotropy and is the gap of spin wave excitation 
that naively could not be bridged by a small deviation 
from the equilibrium point. 

The spin wave dispersion (|12|) has a branch with 
positive imaginary part: u>a{q) — ^(q) at the re- 
gion of wave vectors < \q\ < <?+, where 2Jq\ = 
V'(C22-Cii) 2 +4C 1 2 2 - (C 22 - Cix) and the index A 
stands for advanced. The spin-wave field <p a (t,q) grows 
exponentially with the increment v(q) = Iiau>A{q) reach- 
ing its maximum at \q\ — q m where in the case b) 
2 J Qm — Cn ~ C22 and in the case not b) q m = 0. 
The spin wave distribution function in momentum space: 
N^(t) = <p a (t,q)4> a (t,q), grows accordingly. The spin 
wave energy Eq.(|SJ) is the sum over all modes: £ sw {t) — 
Ylff^sw{t,Q), where using Eq.© we find the energy of 
one spin wave mode: £ sw (t,q) — —KN^(t). 

Growing instability disturbs the distribution of spin 
waves in the momentum space and makes it different 
from the thermal equilibrium distribution in momentum 
space A^2 = T/u>o(q), where the real and positive u->o(k) 
is the spin wave dispersion in the stable state before the 
pulse has arrived. After the pulse the spin wave insta- 
bility is developing and it in turns gives rise to the spin 
wave scattering in the direction that restore the thermal 
distribution and both these processes could be described 
by the Boltzmann kinetic equation: 

dtNf=2v($Ne-I{Nd, (13) 

where v(q) = hnu>A(q) and the relaxation rate is given by 
the Boltzmann collision integral I(N^). In our situation 
obviously the rate of instability and the collision integral 
in the Boltzmann equation are of the same magnitude. 
The leading contribution to the collision integral comes 
from the two magnon scattering with one magnon having 
large momentum p in the incoming and outgoing states 
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FIG. 1: A snapshot of chaotic time dependent domains in 
the developed stage of domain instability in the easy axis: 
E(m) = — X z m 2 , cluster of 96x96 exchange coupled classical 
spins. Non-linearity M 2 (f,t) = 1 caps the growing amplitude 
of the spin waves into a domain picture. Shades of gray repre- 
sent z-component magnetization with black and white being 
spin up and down. 

and second long wavelength magnon with the momentum 
q in the instability domain of phase space: 

J(JV r ) = J W(q,q')(N(q) - N(q))^-, (14) 

where the scattering probability does not depend on mo- 
menta: W(q, q') = W, and it is is small at low tempera- 
tures W ~ (T/J) 5 / 2 . It could be proven that the prob- 
ability of absorption of long wavelength magnon with 
purely imaginary dispersion u>a(q) = iv{q) by a short 
wavelength magnon as well as the probability of emission 
of such a magnon in the process of two short wavelength 
magnons collision are both zero. The probability of two 
long wavelength magnon scattering is proportional to the 
small factor (C/J) 2 and could be neglected. 

Solution of the Boltzmann kinetic Eq . yields: 

Nrft) = exp {(2u{q) - X )t) (N 0q - + $(*)) , (15) 

where $(<) = iV (exp(xt) - 1), x = Wfl is the re- 
laxation rate, £1 = J q<q+ d d q/(2ir) d is the volume of 
instability region in the momentum space and Nq = 
J q<q+ N 0q -d d q/ (2ir) d /£l is the average number of would 
be instable spin waves in the system per one mode before 
the pulse arrives. Solution of Ea. ((r3)) has an asymptote: 
Nq{t) = Noexp(2i/(q)t), as t — ► oo that does not depend 
on the initial occupation number for this mode. The spin 
wave energy Eq.© with 4> a (t) from Eq.© reads: 

£ sw {t) = -V K ( N&)44s, (16) 

Jq<q+ \ Z7T ) 



where V is the volume of the particle and k is the twist of 
precession curve. For large t the spin wave energy could 
be estimated by the steepest descent method. 

Small size of a particle suppresses the domain instabil- 
ity growth. The maximum wavevector of the instability 
region could be shorter then the smallest size-quantized 
wavevector in finite particle. The critical size L c when 
domain instability could develop depends on the shape of 
the particle and is proportional to the domain wall width 
L c ~ y/ J/X, where A is the typical anisotropy energy. We 
studied numerically the case of of 2D square particle with 
periodical boundary condition (torus) with the easy-axis 
perpendicular anisotropy: E(m) = — Arn 2 . For initial 
state n z — 0, we find the critical size L c sw 7^/jJX. 

In conclusion we have found the microscopic equations 
that govern spin wave instability in the non-equilibrium 
precessional state of a ferromagnetic system. In the de- 
veloped stages of instability spin- waves saturate to a time 
dependent chaotic domain pattern. This multidomain 
state Fig.l lives until the intrinsic relaxation of the fer- 
romagnet eventually establishes a single domain ground 
state, randomly if the ground state is degenerate. This 
helps to understand stochasticity of precessional switch- 
ing observed in Ref.[j|. 
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